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ABSTRACT 

One of the active areas of computational aeroacoustics is the application of the 
Kirchhoff formulas to the problems of the rotating machinery noise prediction. The 
original Kirchhoff formula was derived for a stationary surface. In 1988, Farassat and 
Myers derived a Kirchhoff Formula obtained originally by Morgans using modern 
mathematics. These authors gave a formula particularly useful for applications in 
aeroacoustics. This formula is for a surface moving at subsonic speed. Later in 1995 
these authors derived the Kirchhoff formula for a supersonically moving surface. This 
technical memorandum presents the viewgraphs of a day long workshop by the author 
on the derivation of the Kirchhoff formulas. All necessary background mathematics 
such as differential geometry and multidimensional generalized function theory are 
discussed in these viewgraphs. Abstraction is kept at minimum level here. These 
viewgraphs are also suitable for understanding the derivation and obtaining the 
solutions of the Ffowcs Williams-Hawkings equation. In the first part of this 
memorandum, some introductory remarks are made on generalized functions, the 
derivation of the Kirchhoff formulas and the development and validation of Kirchhoff 
codes. Separate lists of references by Lyrintzis, Long, Strawn and their co-workers are 
given in this memorandum. This publication is aimed at graduate students, physicists 
and engineers who are in need of the understanding and applications of the Kirchhoff 
formulas in acoustics and electromagnetics. 

INTRODUCTION 

When Ffowcs Williams and Hawkings published their now famous paper on the noise 
from moving surfaces in 1969 [1], they used a level of mathematical sophistication 
unfamiliar to engineers who would later be the main users of this work. Advanced 
generalized function theory and differential geometry were employed by these authors 
to derive the Ffowcs Williams-Hawkings (FW-H) equation and to obtain some important 
qualitative results in this paper. The subject of generalized functions is very abstract, 
particularly as described in books written by mathematicians. The level of differential 



geometry needed in acoustics is, however, basic and at the level essentially fully 
developed by the end of the nineteenth century. Both of these subjects are not 
emphasized in engineering education. It is possible to teach advanced generalized 
function theory to engineers if some of the abstractions are left out initially. One needs 
to learn how to work with multidimensional Dirac delta functions and their derivatives 
concentrated on moving surfaces, i.e. with support on moving surfaces. This goal can 
be achieved. 

This technical memorandum is on the derivation of the Kirchhoff formulas for moving 
surfaces. The main part of this memorandum is the copies of the viewgraphs based on 
lectures delivered by the author in the Workshop on Kirchhoff Formulas for Moving 
Surfaces at NASA Langley Research Center on February 15, 1995 (see Appendix). 
Attempt was made to present all the mathematical machinery needed in the derivation 
of Kirchhoff formulas. One of the publications of the author [2], NASA TP-3428 (May 
1994), should also be consulted, if needed, to fill in some details. The author and M. K. 
Myers have published two papers on the derivation of Kirchhoff formula for moving 
surfaces [3, 4] which should be easily comprehended by the readers reading the 
material in the Appendix. 

Below we briefly introduce the concept of Generalized Functions. Then we discuss 
the derivation of the subsonic and supersonic Kirchhoff formulas. Finally we make 
some remarks on the development and validation of codes based on the Kirchhoff 
formulas. 


GENERALIZED FUNCTIONS 

Our main reference for this section is NASA TP-3428 [2]. To derive the Kirchhoff 
formulas for moving surfaces, we need to learn how to manipulate multidimensional 
Dirac delta functions and their derivatives. Some knowledge of differential geometry 
and tensor analysis is also essential. In addition to [2], we give some other useful 
references on generalized functions as well as on differential geometry and tensor 
analysis in this paper [5-13]. To learn about generalized functions, we need a change of 
paradigm in the way we look at ordinary functions. Ordinary functions are locally 
(Lebesgue) integrable functions, i.e., functions that have a finite integral over any finite 
interval. This change of paradigm is actually very familiar in mathematics. For 
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example, learning about fractions, negative numbers and complex numbers involves a 
change of paradigm although we are not told that the change is occurring. 

How do we think of an ordinary function /(*)? We think of this function as a table of 
ordered pairs (x, /(*)). A graph of a function is a plot of this table. In generalized 

function theory, we need to work with mathematical objects such as the Dirac delta 
•function" 8(x) with the sifting property 

r ca Hx)S( X )dx = m (i) 

It can be shown that no ordinary function has this property. The Dirac delta function is 
an example of a generalized function which is not an ordinary function. To include <5(x) 
and other such useful but strange objects in mathematics, we change our method of 
thinking about functions as follows. Suppose we take a space of functions D which will 
be called test function space. We will be more specific about D below. Now given an 
ordinary function /(*), let us define the functional 

F[<p] = JT.o/0<k. <p e D. (2) 

If we take the space D large enough, then there is a possibility that the table of 
functional values F[<p ] where <p e D can identify /(*). This is actually true if we take 

the space D as the space of all c°° functions which are identically zero beyond a 
bounded interval, i.e., with compact support. Therefore, the new paradigm of viewing a 
function is: think of the function f(x) in terms of the table |F[0], Q e £>}. We can 

show that this table includes an uncountable number of elements. 

Next, one shows that the functional F[<t>] given by eq. (2) is linear and continuous for 
an ordinary function f(x) [2, 7-9]. We ask whether all continuous linear functionals are 
produced by ordinary functions from eq. (2). The answer is no. For example, the 
functional 


5[0] = 0(0) 0 e D 


( 3 ) 
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is linear and continuous. Therefore, the class of linear and continuous functionals is 
larger than the class generated by ordinary functions through eq. (2). Now, using our 
new paradigm of thinking of a function as a table generated by the functional rule we 
say: 

a generalized function is identified by the table produced using a continuous linear 
functional on space D. 

By an abuse of terminology, we say that: 

generalized functions are continuous linear functionals on space D. 

By this definition the functional in eq. (3) is (represents) the Dirac delta function! Note 
that each continuous linear functional on space D produces (represents, identifies, 
gives) one generalized function. Ordinary functions then become a subset of 
generalized functions called regular generalized functions. Other functions are called 
singular generalized functions. 

Next the operations on ordinary functions are extended to all generalized functions in 
such a way that they are equivalent to the old definitions when applied to ordinary 
functions. To do this, one should write the operation in the language of functionals on 
space D. For example, the derivative of generalized function F[0] is defined by 

F'[0] = -F[0'] (4) 

In this way, many operations on ordinary functions can be extended to generalized 
functions [2, 5-9]. 

Finally, we mention here that the space of generalized functions on D is called £>'. 
For any singular generalized function F[0], we use eq. (2) with a symbolic function 
f(x) under the integral sign. Here the integral does not represent an ordinary integral 
but stands for the rule specified by F[0]. For example, S(x) is a symbolic function 
which is interpreted as follows. Interpret j S(x)0(x)dx as <5[0] = 0(0) for all <f> e D, 

i.e., in our new way of looking at functions as a table of functional values on space D 
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5 (x) = {*(0), M®] • 


(5) 


Of utmost importance to us are delta functions and their derivatives with support on a 
surface f = 0. Here / = /(3c) or / = f(x, t ). We give the following two results [2] 
assuming that |V/| = 1 on / = 0, which is always possible: 

j <p(x)8(f)dx = jy =0 <t> dS (6) 





(7) 


where Hf is the local mean curvature of the surface / = 0 with dS the element of the 
surface area. Also if the function f(x) has a discontinuity across a surface #(3c) = 0 
with the jump defined as 


A/ = f{g = 0 + )-f(g = 0_), (8) 

then 

V/ = V/ + A / Vg 5(g) (9) 

where V/ is the generalized gradient of /(3c) (see [2]). Finally, we mention here that 
the Green’s function method is valid for finding solutions of differential equations with 
discontinuities (weak solutions) provided that all derivatives in the differential equation 
are viewed as generalized derivatives. 


THE KIRCHHOFF FORMULAS FOR MOVING SURFACES 

Assume that /(3c, r) = 0 is the moving Kirchhoff surface defined such that |V/| = 1 
on this surface. Let <p satisfy the wave equation in the exterior Q of / = 0, i.e., 

□ 2 0 = 0 3c e O (10) 

Extend <p to the entire unbounded space as follows, calling the extended function ij> 
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( 11 ) 


0 


x g £2 
x € Q 


The governing equation for deriving the Kirchhoff formula for moving surfaces is then 
found by applying the generalized wave operator (D’Alembertian) to 0 to get [2-4]: 

□ 2 « = "(*« + 7 *»*)«(/) -~ c J t [M„*S(f)\ -V. [**«(/)] (12) 

where M n = v n / c is the local normal Mach number on / = 0, <p n = d<f> / Bn and 
<f> t = d<j> / dt . 

We can now apply the Green’s function method for the wave operator in the 
unbounded space to eq. (11) to find the Kirchhoff formula for subsonically moving 
surfaces [3]. The formula involves a Doppler singularity making it inappropriate for a 
supersonically moving surface. For supersonic surfaces, we derive the Kirchhoff 
formula for an open surface (e.g. a panel). The reason is that the Kirchhoff surface is 
usually divided into panels and the formula is applied individually to each panel. The 
subsonic formula, applies to both open and closed surfaces. However, the supersonic 
formula differs for open and closed surfaces. If the formula for an open surface is 
known, obtaining the formula for a closed surface is trivial. 

The governing equation for deriving the supersonic Kirchhoff formula for a panel is 

□ 2 * - -(<»„ + X -M^ n y{f)s(f) - £ [*■«*"(/)«(/)] -v ■[♦«»(/)«(/)] 

(13) 

where # (/) is the Heaviside function, / is a function such that f > 0 on the panel and 

/ = / = 0 defines the edge of the panel. The derivatives on the right side of eq. (13) 
are brought inside to get three source terms involving H[f^8(f), H[}^8'(f) and 

5(/) 8(f) [4]. The solutions of the wave equation with these kinds of sources are 

given by the author [2]. The Kirchhoff formula for a supersonically moving surfaces 
using the above method was derived and presented by Farassat and Myers [4]. It is a 
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particularly simple and straightforward result and easy to apply. This formula requires 
the mean curvature H F of the surface Z: F(y; x, t) = [ f{y, *)\ ret - We give the formula 

for calculation of H F in the Appendix in terms of the geometric and kinematic 
parameters of the Kirchhoff surface / = 0. 


SOME REMARKS ON DEVELOPMENT AND VALIDATION OF KIRCHHOFF CODES 

The development of a Kirchhoff code requires a good subroutine for retarded time 
calculation if the Kirchhoff surface is rotating. The possibility of multiple emission times 
for a supersonic panel complicates retarded time calculation, particularly for two nearly 
equal emission times. If the Kirchhoff surface is not selected properly for the supersonic 
formula, there is the possibility of a singularity [4]. This singularity can be avoided as 
suggested by Farassat and Myers [4] or by using two different Kirchhoff surfaces for 
different intervals of the observer time. There is a fool-proof test of the Kirchhoff code 
that must not be ignored by code developers. Both of the Kirchhoff formula for moving 
surfaces, as well as that for a stationary surface, are written such that 0 = 0 inside a 

closed surface. Therefore, to test a Kirchhoff code, use a point source inside the closed 
surface and specify <t>, <j> and 0„ analytically on the Kirchhoff surface / = 0. If the 

observer is now put anywhere inside / = 0 and 0^0, then there is a bug in the code. 
One must rule out conceptual misunderstanding of the parameters in the formulation 
first. It is recommended that one should be familiar with the complete details of the 
derivation of the Kirchhoff formulas to avoid conceptual misunderstanding. 

There have been many derivations of the Kirchhoff formula for uniform rectilinear 
motion of the Kirchhoff surface [14, 15]. These formulas do not have the generality of 
Morgans formula derived and rewritten in a new form using modern mathematics by 
Farassat and Myers [3]. Myers and Hausmann [16] were among the first to use the new 
Kirchhoff formula in aeroacoustics. Other researchers include Lyrintzis, Long, Strawn 
and Di Francescantonio [17]. We give separately the publications of Lyrintzis, Long, 
Strawn and their co-workers. 


CONCLUDING REMARKS 

The availability of high resolution aerodynamics and turbulence simulation make the 
Kirchhoff formulas discussed here attractive in aeroacoustics. The mathematics for 
derivation of these formulas have been under development in the last decade and are 
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well within the reach of modern engineers. The final form of the formulas are simple 
and relatively easy to apply. The present paper is written as a guide to understanding 
the mathematical derivation as well as application of these results. 


The viewgraphs in the Appendix give all the necessary mathematical background for 
the derivation of the Kirchhoff formulas. Note that the mathematical part of the 
Appendix is also suitable for understanding the derivation and the solutions of the 
Ffowcs Williams-Hawkings equation. This publication is aimed at graduate students, 
physicists and engineers who are in need of the understanding and applications of the 
Kirchhoff formulas in acoustics and electromagnetism. 

REFERENCES 

1. Ffowcs Williams, J. E., and Hawkings, D. L., "Sound Generated by Turbulence and 
Surfaces in Arbitrary Motion," Philosophical Transactions of the Royal Society, Vol. 
A264, 1969, pp. 321-342. 

2. Farassat, F., "Introduction to Generalized Functions With Applications in 
Aerodynamics and Aeroacoustics," NASA TP-3428, May 1994 (corrected April 1996). 
Download from Langley Technical Report Server, 
ftp://techreports.larc.nasa.goV/pub/techreports/larc/94/tp3428.ps.Z . 

3. Farassat, F., and Myers, M. K., "Extension of Kirchhoffs Formula to Radiation from 
Moving Surfaces," Journal of Sound and Vibration, 123(3), 1988, pp. 451-460. 

4. Farassat, F., and Myers, M. K., "The Kirchhoff Formula for a Supersonically 
Moving Surface," Presented at First Joint CEAS/AIAA Aeroacoustics Conference (16th 
AIAA Conference), June 1995, Munich, Germany. 

5. Lighthill, M. J., Introduction to Fourier Analysis and Generalized Functions, Camb. 
Univ. Press, 1964. 

6. Jones, D. S., The Theory of Generalized Functions, 2nd ed., Camb. Univ. Press, 

1982. 

7. Gel’fand, I. M., and Shilov, G. E., Generalized Functions, Vol. 1 , Properties and 
Operations, Academic Press, 1964. 

8. Kanwal, R. P., Generalized Functions - Theory and Technique, Academic Press, 

1983. 

9. Strichartz, R. S., A Guide to Distribution Theory and Fourier Transforms, CRC 
Press, 1994. 


8 



10. Struik, D. J., Lectures on Classical Differential Geometry, 2nd ed., Dover Books, 
1988. 

1 1 . Kreyszig, E., Differential Geometry, Dover Books, 1 991 . 

12. McConnell, A. J., Applications of Tensor Analysis, Dover Books, 1957. 

13. Aris, R., Vectors, Tensors, and the Basic Equations of Fluid Mechanics, Dover 
Books, 1989. 

14. Lyrintzis, A. S., and George, A. R., "Transonic Blade -Vortex Interactions: The Far 
Field," paper presented at the American Helicopter Society National Specialists’ Meeting 
on Aerodynamics and Aeroacoustics, 25-27 February 1987, Arlington, Texas. 

15. Morino, L., "Boundary Integral Equations in Aerodynamics," Applied Mechanics 
Review, 46, 1 993, pp. 445-466. 

16. Myers, M. K., and Hausmann, J. S., "On the Application of the Kirchhoff Formula 
for Moving Surfaces," Journal of Sound and Vibration, 139(1), 1990, pp. 174-178. 

17. Di Francescantonio, P., "A Supersonically Moving Kirchhoff Surface Method for 
Delocalized High Speed Rotor Noise Prediction," AIAA-96-1708, 2nd CEAS/AIAA 
Aeroacoustics Conference, May 6-8, 1 996, State College, Pa. 

PUBLICATIONS OF A. S. LYRINTZIS AND CO-WORKERS 
Journal/Refereed Papers 

LY1. George, A. R., and Lyrintzis, A. S., "Acoustics of Transonic Blade-Vortex 
Interactions," AIAA Journal, Vol. 26, No. 7, July 1988, pp. 769-776. 

LY2. Lyrintzis, A. S., and George, A. R., "Study of Transonic Blade-Vortex Interaction 
Noise," Noise Control Engineering Journal, Vol. 32, No. 3, May-June 1989, pp. 105-109. 

LY3. Lyrintzis, A. S., and George, A. R., "Far-Field Noise of Transonic Blade-Vortex 
Interactions," American Helicopter Society Journal, Vol. 34, No. 3, July 1989, pp. 30- 
39. 

LY4. Lyrintzis, A. S., and George, A. R., "The Use of Kirchhoff Method in Acoustics," 
AIAA Journal, Vol. 27, No. 10, Oct. 1989, pp. 1451-1453. 

LY5. Lyrintzis, A. S., and Xue, Y., "A Study of the Noise Mechanisms of Transonic 
Blade-Vortex Interactions," AIAA Journal, Vol. 29, No. 10, Oct. 1991, pp. 1562-1572. 

LY6. Lyrintzis, A. S., Wissink, A. M., and Chronopoulos, A. T., "Efficient Iterative 
Methods for the Transonic Small Disturbance Equations," AIAA Journal, Vol. 30, No. 10, 
Oct. 1992, pp. 2556-2558. 


9 



LY7. Xue, Y., and Lyrintzis, A. S., “Transonic Blade-Vortex Interactions: Noise 
Reduction Techniques," AIAA Journal of Aircraft , Vol. 30, No. 3, May-June 1993, pp. 
408-411. 

LY8. Xue, Y., and Lyrintzis, A. S., “Rotating Kirchhoff Formulation for 3-D Transonic 
BVI Noise for a Hovering Rotor, “ AIAA Journal, Vol. 32, No. 7, July 1994, pp. 1350- 
1359. 

LY9. Lyrintzis, A. S., Lee, J., and Xue, Y., “Mechanisms and Directivity of Unsteady 
Transonic Flow Noise," ASME Journal of Fluids Engineering, Vol. 116, No. 3, Sept. 
1994, pp. 649-652. 

LY10. Lyrintzis, A. S., “Review: The Use of Kirchhoffs Method in Computational 
Aeroacoustics," ASME Journal of Fluids Engineering, Vol. 116, No. 4, Dec. 1994, pp. 
665-676. 

LY11. Lyrintzis, A. S., and Mankbadi, R. R., "On the Prediction of the Far-Field Jet 
Noise Using Kirchhoff's Method," AIAA Journal , Vol. 34, No. 2, February 1996, pp. 413- 
416. 

LY12. Lyrintzis, A. S., and Koutsavdis, E. K., “Rotorcraft Impulsive Noise Prediction 
Using a Rotating Kirchhoff Formulation," AIAA Journal of Aircraft Vol. 33, No. 6, Nov.- 
Dee. 1996. 

LY13 Lyrintzis, A. S., Koutsavdis, E. K., and Strawn, R. C., "A Comparison of 
Computational Aeroacoustic Prediction Methods," AHS Journal { in press). 

LY14. Lyrintzis, A. S., and Xue, Y. "Towards a Versatile Kirchhoff’s Method Code," 
AIAA Journal (in press). 

Conference Papers 

LY15. George, A. R., and Lyrintzis, A. S., "Mid-Field and Far-Field Calculations of 
Transonic Blade-Vortex Interactions," AIAA paper 86-1854, AIAA 10th Aeroacoustics 
Conference, Seattle, WA, July 9-11, 1986. 

LY16. Lyrintzis, A. S., and George, A. R., "Transonic Blade-Vortex Interactions: The 
Far-Field," Proceedings of the American Helicopter Society National Specialists’ 
Meeting on Aerodynamics and Aeroacoustics, Arlington, TX, Feb. 25-27, 1987. 

LY17. Lyrintzis, A. S., and George, A. R., "Calculation of Far-Field Noise Using the 
Kirchhoff Method," AIAA Paper 87-2673, AIAA 11th Aeroacoustics Conference, 
Sunnyvale, CA, Oct. 19-21, 1987. 

LY18. Lyrintzis, A. S., and George, A. R., "A Parametric Study of Transonic Blade- 
Vortex Interactions," Proceedings of NOISE-CON 88, Noise Control Design: Methods 
and Practice, edited by J. S. Bolton, Noise Control Foundation, New York, Purdue 
University, West Lafayette, IN, June 20-22, 1988, pp. 155-160. 


10 



LY19. Lyrintzis, A. S., and Xue, Y., "Noise Mechanisms of Transonic Blade-Vortex 
Interactions," Proceedings of the American Helicopter Society 46th Annual Forum, 
Washington, DC, May, 1990. 

LY20. Xue, Y., and Lyrintzis, A. S., "Noise Reduction for Transonic Blade-Vortex 
Interactions," Proceedings of AHS/RAeS International Technical Specialists Meeting on 
Rotorcraft Acoustics and Rotor Fluid Dynamics, Philadelphia, PA, Oct. 1991. 

LY21. Lyrintzis, A. S., Lee, J., and Xue, Y., "Mechanisms and Directivity of Unsteady 
Transonic Flow Noise," presented at the International Symposium on Flow Induced 
Vibration and Noise III, Vol. 3: Flow-Structure and Flow-Sound Interactions, eds 
Farabee, T. M., Paidoussis, M. P., ASME Winter Annual Meeting, Anaheim, CA, Nov. 
1992, pp. 85-113. 

LY22. Xue, Y., and Lyrintzis, A. S., "The Use of Rotating Kirchhoff Formulation for 3-D 
Transonic BVI Far-Field Noise," Proceedings of the American Helicopter Society 49th 
Annual Forum, St. Louis, MO, May 1993. 

LY23. Lyrintzis, A. S., "The Use of Kirchhoff Method in Aeroacoustics," in 
Computational Aero- and Hydro-Acoustics, FED Vol. 147, eds: Mankbadi, R. R., 
Lyrintzis, A. S., Baysal, O., Povinelli, L. A., and Hussaini, M. Y., pp. 53-61; ASME Fluids 
Engineering Conference, Washington, DC, June 1993. 

LY24. Lyrintzis, A. S., Xue, Y., and Kilaras, M. S., "The Use of a Rotating Kirchhoff 
Formulation for the Calculation of High-Speed Impulsive Noise," AIAA paper 94-0463, 
presented at the AIAA 32nd Aerospace Science Meeting, Reno, Nevada, Jan. 1994. 

LY25. Lyrintzis, A. S., Kilaras, M., and Xue, Y., "Transonic 3-D BVI Noise Using a 
Rotating Kirchhoff Formulation for Advancing Rotors," Proceedings of the AHS 50th 
Annual Forum, Vol. I, Washington, DC, May 1994, pp. 115-127; also University of 
Minnesota Supercomputer Institute Research Report, UMSI 94/80, Minneapolis, MN, 
Apr. 1994. 

LY26. Lyrintzis, A. S., "The Use of Kirchhoff’s Method in Rotorcraft Aeroacoustics," 
Paper No. 34, presented at the 75th AGARD Fluid Dynamics Panel Meeting and 
Symposium on Aerodynamics and Aeroacoustics of Rotorcraft, Berlin, Germany, Oct. 
1994; AGARD Conference Proceedings, No. 552, Aug. 1995, pp. 34-1 - 34-16. 

LY27. Lyrintzis, A. S., "A Review of The Uses of Kirchhoff’s Method in Computational 
Aeroacoustics," presented at the 1994 ASME International Mechanical Engineering 
Congress and Exposition, Chicago, IL, Nov. 1994; Acoustic Radiation and Wave 
Propagation, Ed. S. F. Wu, NCA Vol. 17, pp. 73-92. 

LY28. Lyrintzis, A. S., and Mankbadi, R. R., "On the Prediction of the Far-Field Jet 
Noise Using Kirchhoff’s Method," AIAA paper No. 95-0508 presented at the 33rd 
Aerospace Science Meeting, Reno, NV, Jan. 1 995. 

LY29. Strawn, R. C., Biswas, R., and Lyrintzis, A. S., "Computation of Helicopter Blade- 
Vortex Interaction Noise with Kirchhoff Methods" Proceedings of the 51st AHS Annual 


11 



Forum, Fort Worth, TX, May 1995, Vol. I, pp. 495-508; also IMACS Journal of 
Computational Acoustics (in press). 

LY30. Lyrintzis, A., Koutsavdis, E., Berezin, C., Visintainer, J., and Pollack, M., 
"Kirchhoff Acoustic Methodology Validation and Implementation to TiltRotor 
Aeroacoustic Codes (TRAC)," Proceedings of the 2nd AHS International 
Aeromechanics Specialists’ Conference, Vol. I, Bridgeport, CT, Oct. 1995. 

LY31. Lyrintzis, A. S., Koutsavdis, E. K., and Strawn, R. C., "A Comparison of 
Computational Aeroacoustic Prediction Methods," Proceedings of the 2nd AHS 
International Aeromechanics Specialists’ Conference, Vol. I, Bridgeport, CT, Oct. 1995, 
pp. 3-58 - 3-69. 

LY32. Lyrintzis, A. S., "The Use of Kirchhoff s Method in Computational Aeroacoustics," 
presented at the Minisymposium on Issues and Methods in Computational 
Aeroacoustics, 1995 SIAM Annual Meeting, Charlotte, NC, Oct. 1995 [invited]. 

LY33. Wissink, A. M., Lyrintzis, A. S., Oliker, L., Biswas, R., and Strawn, R. C., 
■Efficient Helicopter Aerodynamic and Aeroacoustic Predictions on Parallel Computers," 
AIAA paper No. 96-0153, presented at the AIAA 34th Aerospace Science Meeting, 
Reno, NV, Jan. 1996. 

LY34. Pilon, A., and Lyrintzis, A. S., "On the Development of a Modified Kirchhoff 
Method for Supersonic Jet Aeroacoustics," AIAA paper No. 96-1709, presented at the 
2nd AIAA/CEAS Aeroacoustics Meeting, (17th AIAA Aeroacoustics Meeting) State 
College, PA, May 1996. 

LY35. Lyrintzis, A. S., and Xue, Y. "Towards a Versatile Kirchhoff’s Method Code," 
AIAA paper No. 96-1710, presented at the 2nd AIAA/CEAS Aeroacoustics Meeting, 
(17th AIAA Aeroacoustics Meeting) State College, PA, May 1996. 

LY36. Brentner, K. S., Lyrintzis, A. S., and Koutsavdis, E K., "A Comparison of 
Computational Aeroacoustic Prediction Methods for Transonic Rotor Noise," 
Proceedings of the 52nd AHS Annual Forum, Washington, DC, June 1996 Vol. 2, pp. 
1103-1114. 

PUBLICATIONS OF L. N. LONG AND CO-WORKERS 
Joumal/Refereed Papers 

LOI. V. Ahuja and L. N. Long, "Electromagnetic Scattering using a Parallel Finite 
Volume Runge-Kutta Algorithm," submitted to Journal of Computational. Physics, April, 
1996. 

L02. Y. Ozyoruk and L. N. Long, "Computation of Sound Radiating from Engine 
Inlets, AIAA Journal, Vol. 34, No. 5, pp. 894-901, May, 1996. 


12 



L03. Y. Ozyoruk and L. N. Long, "A New Efficient Algorithm for Computational 
Aeroacoustics on Parallel Processors, Journal of Computational Physics, Vol. 125, No. 
1, pp. 135-149, April, 1996. 

L04. Y. Ozyoruk and L. Long, "Multigrid Acceleration of a High-Resolution 
Aeroacoustics Scheme," to appear AIAA Journal, 1997. 

Conference Papers 

L05. Y. Ozyoruk and L. N. Long, "Progress in Time-Domain Calculations of Ducted 
Fan Noise: Multigrid Acceleration of a High-Resolution CAA Scheme," AIAA Paper 96- 
1771, 17th AIAA Aeroacoustics Conference, State College, PA, May, 1996. 

L06. Y. Ozyoruk and L. N. Long, "A Time-Domain Implementation of Surface 
Acoustic Impedance Boundary Conditions With and Without Flow," AIAA Paper 96- 
1663, 17th AIAA Aeroacoustics Conference, State College, PA, May, 1996. (submitted 
to Journal of Computational Acoustics) 

L07. T. S. Chyczewski and L. N. Long, "Direct Simulations of Supersonic Jet Noise 
on Parallel Computers," AIAA-96-1730, 17th AIAA Aeroacoustics Conference, State 
College, PA, May, 1996. (submitted to Journal of Computational Physics) 

L08. P. J. Morris, L. N. Long, A. Bangalore, T. Chyczewski, D. Lockard, and Y. 
Ozyoruk, "Experiences in the Practical Application of Computational Aeroacoustics," 
ASME Fluids Conference, Invited Paper, April, 1996 

L09. V. Ahuja and L. N. Long, "A Message Passing Scheme for Computational 
Electromagnetics," ACES Conference, Monterey, CA, Mar., 1996. 

LOIO. L. N. Long, "Recent Developments in Computational Aeroacoustics," Invited 
Lecture. 130th Annual Meeting, St. Louis, The Journal of the Acoustical Society of 
America, Vol. 98, No. 5, Pt. 2, Nov., 1995. 

LO1 1. Y. Ozyoruk and L. N. Long, "Computation of Sound Radiation from Engine 
Inlets," CEAS/AIAA Paper 95-063, First Joint CEAS/AIAA Aeroacoustics Conference, 
Munich, Germany, June, 1995. 

L012. V. Ahuja and L. N. Long, "A Finite Volume Algorithm for Maxwell’s Equations on 
Massively Parallel Machines," AIAA-95-1967, AIAA Plasma Dynamics and Lasers 
Conference, San Diego, CA, June, 1995. 

L013. Y. Ozyoruk and L. N. Long, ‘A New Efficient Algorithm for Computational 
Aeroacoustics on Massively Parallel Computers," AIAA-95-1693, Proceedings of the 
AIAA CFD Conference, San Diego, CA, June, 1995. 

L014. V. Ahuja and L. N. Long, "A Message Passing Scheme for Computational 
Electromagnetics," ACES Conference, Monterey, CA, Mar. 20-26, 1995. 


13 



L015. Y. Ozyoruk and L.N. Long, "A Navier-Stokes/Kirchhoff Method for Noise 
Radiation from Ducted Fans," AIAA Paper 94-0462, presented at the 32nd AIAA 
Aerospace Sciences Meeting, Reno, Nevada, Jan., 1994. 

L016. Y. Ozyoruk and L.N. Long, "A Hybrid Scheme for Noise Radiation from Ducted 
Fans," in Proceedings of Noise-Con ’93, Williamsburg, VA, May, 1993. 

PUBLICATIONS OF R. C. STRAWN AND CO-WORKERS 

51. Ahmad, J., Duque, E. P. N., and Strawn, R. C., "Computations of Rotorcraft 
Aeroacoustics with a Navier-Stokes/Kirchhoff Method, "presented at the 22nd European 
Rotorcraft Forum, Brighton, UK, Sept. 16-19, 1996. 

52. Strawn, R. C., Oliker, L., and Biswas, R., "New Computational Methods for the 
Prediction and Analysis of Helicopter Acoustics," AIAA 96-1696, presented at the 2nd 
AIAA/CEAS Aeroacoustics Conference, State College, PA, May 6-8, 1996 (submitted to 
the AIAA Journal of Aircraft) . 

53. Duque, E. P. N., Strawn, R. C., Ahmad, J., and Biswas, R., "An Overset Grid 
Navier-Stokes/Kirchhoff-Surface Method for Rotorcraft Aeroacoustic Predictions," AIAA- 
96-0152, presented at the 34th Aerospace Sciences Meeting, Reno, NV, Jan. 15-18, 
1996. 

54. Wissink, A. M, Lyrintzis, A. S., Strawn, R. C., Oliker, L., and Biswas, R., "Efficient 
Helicopter Aerodynamic and Aeroacoustic Predictions on Parallel Computers," AIAA-96- 
0153, presented at the 34th Aerospace Sciences Meeting, Reno, NV, Jan. 15-18, 1996. 

55. Lyrintzis, A. S., Koutsavdis, E. K., and Strawn, R. C., ‘A Comparison of 
Computational Aeroacoustic Prediction Methods," presented at the 2nd International 
Aeromechanics Specialists’ Conference, Bridgeport, CT, Oct. 11-13 1995, (submitted as 
a Technical Note to the AHS Journal). 

56. Strawn, R. C., Biswas, R., and Lyrintzis, A. S., "Helicopter Noise Predictions Using 
Kirchhoff Methods," presented at the 51st Annual Forum of the American Helicopter 
Society, May 9-11, 1995 (see also Journal of Computational Acoustics, Vol. 4, No. 3, 
Sept. 1996). 

57. Strawn, R. C., and Biswas, R., "Numerical Simulation of Helicopter Aerodynamics 
and Acoustics," presented at the 6th International Congress on Computational and 
Applied Mathematics, Brussels, Belgium, July 25-29, 1994. (see also Journal of 
Computational and Applied Mathematics, Vol. 66, No. 1-2, pp. 471-483 Jan. 1996) 

58. Strawn, R. C., and Biswas, R., "Computation of Helicopter Rotor Acoustics in 
Forward Flight," presented at the 19th Army Science Conference, Orlando, FL, June 11- 
14, 1994. (see also Journal of the American Helicopter Society, Vol. 4, No. 3, July, 
1995, pp. 66-72.) 


14 



S9. Strawn, R. C., Biswas, R., and Garceau, M., "Unstructured Adaptive Mesh 
Computations of Rotorcraft High-Speed Impulsive Noise," AIAA-93-4359, presented at 
the AIAA 15th Aeroacoustics Conference, Long Beach, CA, Oct. 1993. (see also 
Journal of Aircraft, Vol. 32, No. 4, July-August, 1995, pp. 754-760) 


15 



APPENDIX 


Workshop Viewgraphs 


16 



The Kirchhoff Formulas for Moving 
Surfaces in Aeroacoustics — The 



F. Farassat, Aeroacoustics Branch 
Fluid Mechanics and Acoustics Division 
NASA Langley Research Center 
Hampton, Virginia 




© 

JQ 

E 

© 

O. 

© 

CO 


18 


F. Farassat, Aeroacoustics Branch 
Fluid Mechanics and Acoustics Division 
NASA Langley Research Center 
Hampton, Virginia 



Classical Kirchhoff Formula (1882) 
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Classical Kirchhoff Formula (Cont’d) 
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Why are Kirchhoff Formulas Important in Acoustics 
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A Familiar Example of Thinking About Functions by 

New Model 
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Elementary Generalized Function Theory 
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Definition of Generalized Functions (Cont’d) 
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Some Operations on Generalized Functions 
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Some Operations on Generalized Functions (Cont’d) 
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decreasing test functions 5. For our purpose, the most important operation on 
generalized functions is differentiation. 




Differentiation of Generalized Functions 

All test functions are in D. 

• f(x ) ordinary function, differentiable, F[< |>] = J fydx, we must identify F'[<|)] 
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Some Important Results of Generalized Function 
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Green’s Function of a 2nd Order Linear O.D.E. 
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Green’s Function of a 2nd Order Linear O.D.E. 
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Generalized Functions in Multidimensions 
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Some Results From Differential Geometry (Cont’d) 
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Some Results From Differential Geometry (Cont’d) 

The principal directions can be found from solving the quadratic equation: 
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Some Results From Differential Geometry (Cont’d) 
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Some Results From Differential Geometry (Cont’d) 
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Integration of Product of Delta Functions (Cont’d) 
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Illustration of Manipulation of Generalized Functions 

(Cont’d) 

In deriving conservation laws in differential form from finite volumes involving 
discontinuous functions, whenever the divergence theorem is used to convert 
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Things to Know About Green’s Function 
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Use of Green’s Functions for Discontinuous Solutions 
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The Two Forms of the Solution of Wave Equation 

(Volume Sources) 
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The Two Forms of the Solution of Wave Equation 

(Volume Sources) (Cont’d) 
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The Governing Wave Equation for Deriving 
Kirchhoff Formulas (Cont’d) 
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Derivation of the Classical Kirchhoff Formula (Cont’d) 
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Derivation of the Subsonic Kirchhoff Formula 
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Derivation of the Subsonic Kirchhoff Formula (Cont’d) 
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Derivation of the Subsonic Kirchhoff Formula (Cont’d) 
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Derivation of the Subsonic Kirchhoff Formula (Cont’d) 
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Derivation of the Subsonic Kirchhoff Formula (Cont’d) 

Similar procedure for / 2 gives 
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Derivation of the Subsonic Kirchhoff Formula (Cont’d) 
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A Simple Trick in Preparation for Supersonic 

Kirchhoff Formula 
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A Simple Trick in Preparation for Supersonic 
Kirchhoff Formula (Cont’d) 
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Derivation of the Supersonic Kirchhoff Formula 
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Derivation of the Supersonic Kirchhoff Formula 

(Cont’d) 

Now putting the solutions for , <j>2 ? and together in (j), we get the 
supersonic Kirchhoff formula 
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Selection of the Kirchhoff Surface for Supersonic 

Kirchhoff Formula 
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References for Generalized Functions 



© 

-Q 

E 

© 

Q. 

a> 

co 


F. Farassat, Aeroacoustics Branch 
Fluid Mechanics and Acoustics Division 
NASA Langley Research Center 
Hampton, Virginia 




E 

<D 

Ql 

<D 

CO 


F. Farassat, Aeroacoustics Branch 
Fluid Mechanics and Acoustics Division 
NASA Langley Research Center 
Hampton, Virginia 



REPORT DOCUMENTATION PAGE 

Form Approved 

OMBNo. 0704-0188 

Pubic iipotling buidan for tots ocaeceon of aitemabon t» eabmated to average i hour per response, includfog the firm tor iwaewring inatnietions, aaarahingeiutirig data eou'cet, 
gatoaring and maintaining 4 m data needed, and complefing and reviewing the collection of information. Send oommenti regandng Bus burden ertmato or my other aspect of tha 
oofiacbon of ^formation. tndudmQ tugge aborts tor rxfaoog the burden, to Washington 1 leadguartare Services. Directorate for Information Operations and Reports, 1215 Jeftareon Davis 
Highway, Suite 1204. Arlington. VA 22202-4302, and to toe Office of Management and Budget Paperwork Reduction Project (0704-0188). Washington, DC 20603. 

1. AGENCY USE ONLY (Lmvw blank) 

2. REPORT DATE 

September 1996 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

The Kirchhoff Formulas for Moving Surfaces in Aeroacoustics-The 
Subsonic and Supersonic Cases 

5. FUNDING NUMBERS 

WU 538-03-12-01 

6. AUTHOR(S) 

F. Farassat 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 


0. SPONSORING / MONITORING AGENCY NAJIE(S) AND ADORESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING / MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 110285 


11. SUPPLEMENTARY NOTES 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 71 

Availability: NASACASI, (301) 621-0390 


12b. DISTRIBUTION COOE 


13. ABSTRACT (Maximum 200 word*) 

One of the active areas of computational aeroacoustics is the application of the Kirchhoff formulas to the problems of the 
rotating machinery noise predictions. The original Kirchhoff formula was derived for a stationary surface. In 1988, Farassat 
and Myers derived a Kirchhoff Formula obtained originally by Morgans using modem mathematics. These authors gave a 
formula particularly useful for applications in aeroacoustics. This formula is for a surface moving at subsonic speed. Later in 
1995 these authors derived the Kirchhoff formula for a supersonically moving surface. This technical memorandum presents 
the viewgraphs of a day long workshop by the author on the derivation of the Kirchhoff formulas. All necessary background 
mathematics such as differential geometry and multidimensional generalized function theory are discussed in these 
viewgraphs. Abstraction is kept at minimum level here. These viewgraphs are also suitable for understanding the derivation 
and obtaining the solutions of the Ffowcs Williams-Hawkings equation. In the first part of this memorandum, some 
introductory remarks are made on generalized functions, the derivation of the Kirchhoff formulas and the development and 
validation of Kirchhoff codes. Separate lists of references by Lyrintzis, Long, Strawn and their co-workers are given in this 
memorandum. This publication is aimed at graduate students, physicists and engineers who are in need of the 
understanding and applications of the Kirchhoff formulas in acoustics and electromagnetics. 


14. SUBJECT TERMS 

Aeroacoustics, Kirchhoff Formulas, Noise Prediction, Generalized Functions 


17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

Unclassified 

Unclassified 

Unclassified 


15. NUMBER OF PAGES 
112 


16. PRICE CODE 

A06 


20. LIMITATION OF ABSTRACT 

Unlimited 


Standard Form 296 {Rev. 2-89) 
Prescribed by ANSI Sid 239 1 8 


NSN 7540-01-280-5500 














